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Abstract 

We present a general analysis on charmless B-meson decays B —* tttt and 
ttK. It is noticed that the final state interactions and inelastic rescattering 
effects must be significant in order to understand the consistency of the current 
data. By using general isospin decompositions, the isospin amplitudes and 
the corresponding strong phases could be extracted from a global x 2 fit of the 
experimental data. We emphasize that in general there are two, rather than 
one, relative strong phases in the decomposition. In the assumption of two 
equalstrong phases as considered in the literature, the current data , especially 
the ones concerning B — ► ^i^Wdecays, will imply a large isospin amplitude 
l a 3/2l ^ ^O, which is larger by a factor of five than the one from the naive 
factorization estimation. When two different strong phases are considered, all 
the isospin amplitudes can become, within the lc-level, comparable with the 
theoretical values. We also show that the difference between the two strong 
phases cannot be too large and will be restricted by the most recent upper 
bound of B — > 7r°7r° decay. In any case, the strong phases are found to be large 
and the branching ratio of B — > 7r°7r° is likely to be enhanced by an order of 
magnitude in comparison with the one obtained from the naive factorization 
approach. Direct CP violations in all decay modes are also calculated and 
found to be close to the sensitivity of the present experiments. 

PACS 11.30.Er, 12.60.Fr 
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I. INTRODUCTION 



Recently the CLEO collaboration has reported measurements on the branching ratios 
of rare hadronic B decays B — > 7111,71 K [1|,[|. The data have attracted great interest from 
both theorists and experimentalists. The study of these channels will provide us impor- 
tant insights on understanding the effects of electroweak penguins (EWP) in B-system 
@, and the final state interactions(FSI) ||||, as well as extracting the weak CKM phase 
7 = arg(y u dV* b /V c dV* b ) P-pyiOH . It may also open a window for probing new physics [[LJ] . 

From the current data, it is noticed that the branching ratio for B — > 7r + 7r~ is relatively 
small, Br(B — > 7r + 7r~) ~ 4 (in units of 10~ 6 ). The decays for B — > tt + K~ , tt~K° have almost 
equal decay rates, i.e., BriB — > tt + K~) ~ BriB — > 7r~K°) ~ 17. While an unexpectedly 
large branching ratio for B — > tt°K° decay was also observed, BriB — > 7t°K°) ~ 14. 
These measurements seem in conflict with the calculations based on the naive factorization 
hypotheses. For the first three decays, it was pointed out that the factorization approach 
may still be valid if one takes the weak phase 7 to be greater than 90° ||. With such a large 
7, i.e., cos 7 < 0, the interference between tree and penguin diagrams has opposite sign in 
B — > mi and B — > ttK decays. Thus the negative cos 7 will suppress the decay rate for 
B — > 7T7T and enhance that for B — > it + K~ . As a consequence, the almost equal decay rates 
for B — > tt + K~ and B — > tt~K° decay modes indicate the dominance of strong penguin. 
However, the large rate for B — > 7t°K° is not easily explained. Most recent analysis showed 
that a large final state interaction (FSI) phase would be helpful to enhance the branching 
ratio for the B —>■ tc°K° decay [|j], but only considering the elastic- rescatterings remains 
insufficient to obtain the large central value of the data. 

To understand the measured data, besides some model-dependent calculations, a model- 
independent approach using a single relative strong phase has also been proposed for the 
study of B — ► 7T7T, ixK decays f§|5||§ . Since the ordinary factorization approach suffers from 
the uncertainties due to hadronic matrix elements, such as the meson decay constants, B- 
meson form factors and so-called effective color number N c . The mo del- independent analyses 
may be more useful as more data become available. The approach based on the isospin SU(2) 
and approximate flavor SU (3) [0 symmetries of the strong interactions has been proposed 
to constrain and extract the weak phase 7 [0]. It has been noticed that the ratios 

between the CP-averaged decay rates, such as R = Br[B — > n ± K°)/Br(B — > ^K^) may 
provide us important information on the weak phase 7. Most recently, it has been shown 
that the weak phase 7 may be determined through three ratios among CP-averaged decay 
rates of B 7r + 7r~, 7r + 7r°, n~K + and tt~K° decays ||I0 |. Where two solutions were obtained 
at the la level, one with positive cos S and negative cos 7, i.e., relative small strong phase 
S and large weak phase 7, and another with negative cos 5 and positive cos 7. The latter 
with positive cos 7 seems to be favored by solutions obtained from other constraints in the 
standard model but appears not to be as favorable as the one with negative cos 7 studies of 
all the existing charmless decays are taken in account. However, there is still no complete 
analysis in the literature. 

In this paper, we shall give a general analysis for all 7 decay modes of B — ► tttt and 7r.fr. 
For that purpose, we will start from a general model independent parameterization for all the 
decay amplitudes by considering both the isospin and simple diagrammatic decompositions. 
We will show that there are in general 15 independent variables. By assuming the SU(3) 
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relations which appropriately account for SU(3) symmetry breaking effects, it allows us 
to reduce the 15 independent variables into 9. They consist of 6 isospin amplitudes and 2 
relative strong phase as well as one weak phase 7. Note that once the relative strong phase is 
zero, B — > tt~K° only receives contributions from penguin- type diagram and the amplitudes 
of the isospin I = 1/2 and 1 = 3/2 amplitudes from tree- type diagrams cancel each other, 
which provides an additional constraint ||. In fact, one of the isospin amplitudes becomes 
almost irrelevant due to the suppression factor of the CKM mixing element. With such a 
consideration, there are only 8 relevant unknown quantities with 5 isospin amplitudes and 
3 phases. We show that the current 6 measured decay rates allow us to extract 6 unknown 
quantities as functions of two variable. The upper bound of the decay rate B — > 7r°7r° also 
provides a bound for the difference between the two strong phases. Once taking the numerical 
value of the weak phase 7 to be the one obtained from other constraints in the standard 
model and fixing one of the strong phases, all the other parameters can be determined. 
With these determined parameters, we are then able to predict the branching ratio of the 
B — > 7r°7r° decay mode which is yet unmeasured due to the difficulty of its identification by 
the current detector. In addition we also present predictions for direct CP violations in all 7 
decay channels of B — ► titt, rrK. In our numerical fitting, we have adopted the x 2 analysis 
for the CLEO data in order to have a systematic treatment on the experimental errors. 

In general, according to the Watson theorem, there are two independent relative strong 
phases associating with the isospin amplitudes. They are often assumed to be equal in the 
literatures |4]||||. In this work, we shall make a more general analysis with two relative 
strong phases. It is shown that the equal phase assumption will result in large enhancement 
of isospin amplitude dy 2 which will be 5 times larger than the one calculated from the 
factorization approach. The value of the strong phase is found to be 5 ~ ±95°. These large 
values may imply large inelastic FSI or indicate the possible new physics effects. However, 
if the two strong phases are different, the value of (2*3 /2 Cclll be lower and is comparable with 
the usual factorization calculations. 

It is remarkable to observe that within la all the 6 decay rates can be consistently fitted 
for a large range of the weak phase 0° < 7 < 180° for the above two cases. It is also of 
interest to note that one of isospin amplitudes and the strong phases have a weak dependence 
on the weak phase 7. Three isospin amplitudes show a moderate dependence on the weak 
phase 7. Only one isospin amplitude is sensitive to 7. In particular, the fitting values for 
the 4 usual isospin amplitudes considered in most of the literature could still be comparable 
with the ones obtained by using naive factorization approach. The resulting large strong 
phases may be regarded as a strong indication of large FSI in B — > nn, irK decays. 



II. THE GENERAL FRAMEWORK 



We begin with writing the decay amplitude of B — > ttti, ttK in the following general form: 

j^tttt(ttK) ^d(s) j^tttt(ttK) _|_ yd(s) j^tttt{itK) /-Q 

where = V ub V* d ^ and X^ = V cb V* d ^ are the products of CKM matrix elements. The 

term proportional to Xf^ = V tb V* d ^ has been absorbed into the above two terms by using 
the unitarity relation, V ub V* d{s) + V cb V* d{s) + V tb V t d{ s ) = 0. 
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we also find it useful to adopt the isospin decomposition for the decay amplitudes 



^v^y^o-v^ + y^re* (2) 

AXo = ^o^ - ^are^ (3) 

= (4) 

Ki K - = {i^ 112 + yf «3/ c 2 ^ 3/2 (5) 

Koko = ^l</2^ ~ \Jl<f 2 e^ (6) 

K'ok- = -^<h%^ - ±=A u f_ R0 (7) 



with 



K-ko = f^if/^ - ^%e^ (8) 

Where a"' c and b u { c are the isospin amplitudes and 5i and 8y 2 are the strong phases due 
to final state interactions. In some literatures the strong phase of isospin amplitude by C 2 is 
assumed to be equal to the one of a u ' c for simplicity (|,|8]fJ. However, in the most general 
case, these strong phases are not necessarily the same, since they arise from the effective 
Hamiltonian with different isospin. The subscript 1 = 0, 2, 1/2, 3/2 denote the isospins of 
the amplitudes. The advantage of the isospin decomposition allows one to use SU (3) rela- 
tions including leading order SU (3) breaking effects. In other words, the isospin amplitudes 
are assumed to satisfy the following relations: 

U,C ( p I p \ U,C U,C ( p IP \ u,c 

«o - (/*•/ /*>i/2> a 2 > ~ {U/fK)a 3 ' /2 

$0 - $1/2, $2 ^ $3/2 (9) 

where f n and fx are the it, K meson decay constants with j^j f% ~ 0.8. For convenience, 
we define two phase differences as follows: 

$ = $3/2 ~ $1/2 

$' = $[/2-$l/2 (10) 



Practically, the decay amplitudes are evaluated by calculating various Feynman diagrams. 
In order to see how those isospin amplitudes receive contributions from diagrams, we also 
present a simple diagrammatic decomposition. The diagrams can in general be classified into 
6 types denoted by T(tree), C(color suppressed tree) P(QCD penguin), Pg^electroweak 
penguin) and P^wi c °l°r suppressed electroweak penguin) [|l"5|j : 
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A v - v+ =T + P+lpgy (11) 

A n0n0 = i=(-C + P-P EW - ^P° w ) (12) 
At-tt — TT^ - — ^ — ^bw ~~ Pew) (13) 

A^ K - = r + p' + lp% v (14) 

A^o^o = ~y|( — ^' + -P' — P'ew ~ ^P'ew) (15) 

where the primed and unprimed quantities are the amplitudes in B — > ttK and B — ► 7r7r 
decays. They roughly differ by a factor f^/fx — 0.8 when the SU(3) flavor symmetry 
breaking effects are considered. 

Combining the two decompositions, it is straight forward to get the following relations: 

a$e"v a = -L (2T' - C" + 3P' - P' EW + , (16) 

<h%e^ = ^(T' + C' + P' EW + P' E C W ) U ' C . (17) 

From the above equations, one may easily see the relative magnitudes among those 577(2) 
invariant amplitudes. If the inelastic rescattering effects are small, T', C will only contribute 
to the term proportional to A*. Therefore one may expect that T'(C') U ^> T'(C') C . This 
will lead to ay 2 (3/2) a i/2(3/2)- Since a\, 2 receives contributions from QCD penguins, while 
a^, 2 only gets contributions from EWPs, one may conclude that a\, 2 ^> a\, 2 . 

To obtain relations for the isospin amplitude by C 2 , one needs to be careful in adopting the 
diagrammatic decomposition implied by the naive factorization ansatz. This is because the 
resulting relative strong phase is zero in the factorization approach , i.e., 5 = 63/2 — = 
and 5' = 0. As a consequence, 

ArOK- = -^(-T' -C'-P'- P' EW - \P'C W ) (18) 

K-k« = P'~ \p'ew (19) 

The amplitudes with isospin 1=1/2 and 1 = 3/2 from tree-type graphs cancel each other 
in Eq.(|IT5|). Thus the total amplitude only receives contributions from penguin diagrams in 
this case, namely 

K-ko = ffT/l - f^; 2 = {p> - \p' E c w ) u ' c (20) 

If assuming the i-quark dominance in the penguin loops, one finds from eq. (120) that : 



A^-KO — 0r &T/2 — &1/2 + ~y| a 3/2 _ ~y| a 3/2- ( 2 1) 
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which may be assumed for simplicity to be approximately valid after considering final state 
interactions with nonzero strong phases. In the numerical calculations, we have check that 
the amplitude 6"/ 2 is less important due to the strong suppression of the CKM factor (for 
instance, even if taking b\i 2 ~ b\, 2 , the results remain almost unchanged). 

With the above analyses, let us provide an intuitive discussion on how to yield a large 
branching ratio for B — ► tt°K decay by appropriately choosing the isospin amplitudes. Note 
the fact that A* <C A*, one may roughly estimate the ratio between Br(B — * n°K°) and 
BriB — > 7i + K~) by neglecting the terms containing the CKM factor A^ 



Br(B -> n°K°) 
~ Br(B -»• TT+K-) 



3 a l/2 V 3 a 3/2 ei<5 



1 n c -I- 4 / -n c P iS 
3 u l/2 ^ V 3 3/2 e 



(22) 



with 5 = 5 3 / 2 — 5i/2- Once neglecting a^, 2 , the ratio may be simply given by R ~ |, which 
is much smaller than the central value of the data R = 0.84. It indicates that to enhance 
the decay rate of B — > ir°K°, the isospin amplitude ay 2 should not be neglected. Its small 
value may provide a sizable contribution for a large value of 5 > tt/2. This is because in 
this case there exists a constructive interference between a\i 2 and a^, 2 in B — > it K and a 
destructive interference in B — > 7r + ir~. The situation is quite similar to the case for a large 
7 > 7r/2, which is considered to enhance B — > 7rA" and decrease B ^ tttt decay rates. From 
eq.fl22|), it is easily seen that the value of satisfies 



a c 3 /2> ^ + ^ a; /2 ~0.12 a; /2 . (23) 

With the above considerations, there are only 8 relevant quantities: 
a i/2' a 3/2' ^l/2' ^> 7, which should be constrained by six measured decay rates and 
one upper bound. When taking the weak phase 7 as a free parameter, the rest 6 variables 
can be determined from 6 equations of eqs.(0-0). As the errors in the current data remain 
considerable large, one may not take the central values of the data to be too serious. Thus 
by only using the central values of the data to determine the 6 variables may not be good 
enough. To take into account the experimental errors in a systematic way, we shall adopt 
a global x 2 { least square) analysis for the present data. This treatment allows us to obtain 
not only the central value but also the errors for the fitting amplitudes. Our fitting will be 
carried out by using the standard \ 2 analysis program package MINUIT [TH . 



III. RESULTS AND DISCUSSIONS 

In order to compare with the values estimated from the factorization, it is necessary 
to explicitly see how large of the isospin amplitudes dy C 2 3 / 2 and 61/2 from the factorization 
calculations, we present the relevant formulae for the B — > nn, ttK decay amplitudes with 
the assumption of factorization |J[8[] : 



G 

A n+K ~ = i-^=f K F^(m 2 K )(m 2 B - ml) (A>i - X s t (a 4 + a 10 + (a 6 + a 8 )R 4 )) (24) 

G / 1 

A^oko = -i^fKF^(m 2 K )(m 2 B - ml)X s t (a 4 + a 6 R 5 - - (a 10 + a 8 R 5 
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W°K 



A c - 

7T — 



-i^UF BK (ml)(m 2 B - m 2 K ) (a> 2 - A^(a 9 - a 7 )) 
-i^f K F^ n (m 2 K )(m 2 B - m*) (A>i - A t s (a 4 + a 6 i? 4 + (aio + a 8 R 4 )) 
-i^UF BK (ml)(ml - m\) (a> 2 - A^(a 9 - a 7 )) 
= -i-Zf K F**(m 2 K )(m 2 B - m*)A t s (a 4 + a 6 R 5 - ^(a 10 + a 8j R 5 )) (27) 



(25) 



(26) 



where f nj K and F BlT > BK are the decay constants and 5-meson form factors respectively, 
R4 = 2m 2 K /{rrib — m u )(m u + m s ) and R$ = 2m 2 K /{rrib — md){jrid + m s ). In the flavor SU (2) 
limit, one has i? 4 ~ i? 5 ~ 2m 2 K /{m b m s ). 

The expressions of the isospin amplitudes can be rewritten as follows: 



a u 1/2 e iS ^ 



/2 /I 1 

'' I a/ — (0.1 + a 4 + aio + (a 6 + a 8 )i? 4 ) + W -(a 4 + a 6 i? 5 - -(ai + a 8 i? 5 )) 



-J^(a 2 + ^(a 9 - a 7 )X) 



(28) 



^2 [T 1 

-(a 4 + aio + (a 6 + a 8 )i? 4 ) + y-(a 4 + a 6 i? 5 - -(ai + a 8 i? 5 )) 



1^ 



(o 9 - a 7 )X 



/ 



(29) 



a c 3/2 e^ 



y 3^"+^- ~~ y z a I°k° 
IT ( 3 3 

V 3 V 1 + + 2^ 9 ~~ + 2^ 10 + a8 ^ 4 



'Ha- 



r ^-( a io + a 8 -R 4 + (a g - a 7 )X) 



(30) 



(31) 



and 



= r ^a 4 + a 6 i? 5 - i(a 10 + a 8 R 5 )j 



(32) 
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where r = i^f K F B *{rn 2 K ){ml-ml) and X = (f n /f K ) {F BK /F**){m 2 B -m 2 K )/{m% -ml). 
In our numerical estimates, we will take f n = 133 MeV, fx = 158 MeV, F Bn = 0.36 and 
F BK = 0.41. There remains a large uncertainty in strange quark mass m s . For m s = (100 ~ 
200) MeV, we find that the numerical values of those amplitudes are given by 

a u l/2 ~ 818 ~ 846, a u z/2 ~ 709 , 
a c 1/2 ~ -(103 ~ 131), 4 /2 ~ -7 , 

~ -(72 ~ 100) . (33) 

The results from the x 2_ fitting are shown in Fig.|l}|| where the six amplitudes as well 



as their errors at la level are obtained as functions of the weak phase 7 with 5' fixed at 
0,7r/6,7r/3 respectively. The relative magnitudes of the amplitudes are consistent with the 
previous discussions. In our fit, the minimum value of \ 2 is found to be extremly low ( 
typiclly Xmin ~ 0-5 x 10~ 12 ). This means that the x 2- fits are highly consistent and the 6 
amplitudes are actually extracted as the solutions of Eqs.2-8. It can be seen from the figures 
that the 7 dependence of a\j 2 is quite strong and the one of a\i 2 and a^f 2 is relatively weak. 
On the other hand, it may be used to extract the angle 7 once one of those amplitudes can 
be determined or calculated in other independent ways. It is of interest to see that the 7 
dependence of the amplitude b\, 2 and the strong phase 5 is weak, which shows that these 
two quantities are approximately fixed. The possibility of large 5 was also suggested in Ref. 
Q to explain the large branching ratio of the B — > tc°K decay. Recently, perturbative QCD 



calculations have also shown a large strong phase |L7 . 

In Fig.|l| where the phase difference 5' is set to be zero as usual, the x 2_ fitting shows 
that the values of a\i 2 and a^ 2 may be comparable with the ones from the theoretical 
estimations, only when the weak phase 7 is large. Especially for 7 > 2n/3, the fitting values 
could coincide with the ones from factorization except for a^ 2 . For 7 < 27r/3, the two 
amplitudes are smaller than the ones from naive factorization calculations. It appears that 
the factorization approach may become suitable for large weak phase 7. This phenomenon 
was observed by most of the analyses in the literature which neglects the isospin amplitude 
03/2- As a consequence, the resulting large value of 7 seems to be conflict with the one 
obtained from other constraints in the standard model. Before drawing the final conclusion, 
one also notice that in the factorization approach, one yields a zero strong phase 5 = 0° 
which actually contradicts with the general fitting value 5 ~ ±95°. Therefore, the results 
of estimates based on the naive factorization approach should be unreliable, and the isospin 
amplitudes must receive additional large contributions. A large value for the relative strong 
phase 5 ± 95° implies that the final state interactions or inelastic rescattering effects must 
be significant. 

We would like to stress that the most outstanding feature of the x 2 analysis with 5' = 
is that the isospin amplitude is likely to be relatively larger than the one estimated 
from the naive factorization calculations. The fitting central value of a^, 2 may be larger by 
a factor of 7 ~ 9. To explicitly see how the decay rates depend on the isospin amplitude 
a 3/2' we plot in Fig.[| the 6 branching ratios of B — > nn, ttK decays as functions of a\, 2 . It 
can be clearly seen that if 5' = 0, a small value of ~ —7 is not able to reproduce all the 
CLEO data within the la level, especially the data for the channels of 7r°K° and tt°K~ . To 
consistently describe the whole data, we need a relative large value a^, 2 ~ —75 for fitting 
the central value of the data which is about 10% of the largest one 0%. 
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Within the standard model, it seems difficult to enhance the isospin amplitude a\, 2 by 
an order of magnitude even the inelastic FSI is involved, this is because the main inelastic 
channels such as B — > DD(DD s ,rj c K) — > -ktt(K), only contribute to isospin-| part of the 
decay amplitude. In the SM, it is known that the ratio 03/9/^3/2 can be determined without 
the hadronic uncertainties in the flavor SU(3) limit, this is because the ratio only depends 
on the short distance Wilson coefficients ||19|| . Thus a large value of a>%/2/ a 3/2 ma y indicate 
the existence of new physics. While all models beyond the standard model must effectively 
provide large contributions to the electroweak penguins in order to enhance the isospin 
amplitude cty 2 , such models are: SUSY with R parity violation, Z'-model and Z-mediated 
FCNC models, et. al. 

Let us now consider the case that 5' is non-zero, the situation then becomes quite dif- 
ferent. In Fig.0 and Fig.||, it is seen that a^ 2 decreases as the value of 5' increases. When 
5' reaches ir/3 , a\, 2 will be consistent with the value yielded from factorization approach. 
It is also noticed from the figures that a large 5' leads to large values of a\j 2 and 5. The 
enhancement of a^ 2 may be easily understood than the enhancement of a^ 2 due to final 
state interactions. Nevertheless, as will be discussed below, the values of 5' cannot be too 
large due to the constraint of the upper bound of the branching ratio of B — > 7r°7r°. 

When all the isospin amplitudes and strong phases are determined, one is able to predict 
the direct CP asymmetries for all the relevant decay channels. The direct CP asymmetry 
in B — > 7T7T, 7rK decays is defined in the standard way 

T(B f) - V(B f) f , . 

where / denotes the final state mesons. In Fig.|^ we plot several |Acp|'s as functions of 
the weak phase 7 with different value of 5'. When 7 is near 90° and 5' = one has 
\Acp{'k + K~) \ ~ 0.04 which is in a good agreement with the most recent CLEO data Aqp = 
—0.04 ± 0.16 0. At this point, we have a reliable prediction for the direct CP violations in 
the following decay modes. The Acp's with 45° < 7 < 95° read 

\A CP (7T + K~)\ ~ (2.5 ~ 4)% , \A CP {^K°)\ ~ (2.5 ~ 5)% , 
\A CP (n°K-)\ ~ (7.5 ~ 10)% , \A CP (7i-K )\ ~ (5 ~ 6.5)% , 



|A C p(7r+7r-)| ~ (7.5 ~ 12.5)% , |v4 C p(7r 7r°)| ~ (7.5 ~ 12)% . (35) 
for 5' = 0, and 

\A cp (tt + K-)\ ~ (5 ~ 8)% , \A CP (7r K°)\ ~ (8 ~ 10)% , 

\A cp {tt°K-)\ ~ (12 ~ 18)% , \A CP {n-K°)\ ~ (8 ~ 12)% , 

|A C p(7r + 7r-)| ~ (16 ~ 24)% , \A CP (n°n )\ ~ (10 ~ 14)% . (36) 

for 5' = 30°. In spite of the large 5, the smallness of CP asymmetries in B — » 7t~tt° decay is 
due to the absence of the interference between tree and penguin diagrams. 

There remains an unobserved decay mode in B — » titt, ttK decays, which is the decay 
mode B — > 7r 7r ( the CLEO collaboration has already reported the indication of BriB — > 
7r°7r~) ~ 5.6 ). As all the relevant isospin amplitudes have been determined as functions of 
7, it allows us to predict B — > 7r°7r° as a function of 7. It is interesting to note that our 
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X 2 -analysis shows that the resulting branching ratio BriB — > 7r°7r°) is almost independent 
of the weak phase 7. Its value at 5' = is close to the one of B — > 7r + 7r~ decay 

Br(B -> 7T°7r ) ~ 4.6 x 1(T 6 . (37) 

with 5' increasing, the branching ratio becomes larger and can reach ~ 7(10) x 10~ 6 when 
5' = 30° (60°). Such a large branching ratio is about an order of magnitude larger than 
the prediction based on factorization calculations. Most recently, the CLEO collaboration 
reported an upper bound of BriB — > 7T 7r°) < 5.7 x 10 ( ^ 6) f2J, this will impose a strong 
constraint on the value of 5' (see Fig. 6). It is seen that to be consistent with the data at 
the la level, the upper bound of the branching ratio BriB — > 7r°7r°) limits 5' to be less than 
~ 50°. 

IV. CONCLUSIONS 

In summary, we have made a general less model-dependent investigation on the charmless 
B-meson decays by using the y 2 analysis based on the most recent CLEO data. We have 
used the most general isospin decomposition with two independent strong phases. All the 
isospin amplitudes in rare hadronic B decays B — > nn, nK can be determined as functions 
of the weak phase 7 and one strong phase 5'. The effects of two equal and unequal strong 
phases are studied in detail. It is found that the isospin amplitude b\, 2 and the strong 
phase 8 only slightly depend on the phase 7. An important observation under the equal 
strong phase assumption is the relative large isospin amplitude ay 2 , (a^ 2 ~ —70) where 
the central value is about 5 times greater than the one obtained from the factorization 
calculations. When the two strong phases are not equal, the allowed values of a\, 2 decrease 
as their difference, i.e., 5', increases. For the most general case with two rather than one 
large FSI strong phases, the magnitude of all the isospin amplitudes may be around the 
one estimated from the factorization approach. Nevertheless, one needs to find out the 
mechanism of producing large strong phases. This could directly be tested by measuring 
the branching ratio BriB — > 7r°7r°). 

The direct CP asymmetries Aq P for all the relevant decay channels have also been given 
as functions of 7. The resulting numerical value for AqjF is consistent with the most 
recent data. When taking 7 to be in a reasonable range 7 = 45° ~ 95°, we are led to the 
results given in eqs.(35) and (36), which can be directly tested by experiments in the near 
future. A resulting large branching ratio BriB — > 7r°7r°) which is comparable with the one 
BriB — > 7r + 7r~) will also provide an important and consistent test. 

From the most general analysis presented in this paper, the data appears to strongly 
suggest that final state interactions and inelastic rescattering effects must be significant 
and play an important rule in the charmless B — > ttti, nK decays. Otherwise, our general 
analyses may be interpreted as hinting at the existence of new physics. For a more definite 
conclusion, one needs more precise data. The two B-factories BaBaR and BELLE are 
expected to provide us with more information from the charmless decays. 
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FIG. 1. The isospin amplitude fitted as a function of weak phase 7 from a x 2 ~ analysis of the 
recent CLEO data. The vertical bars indicate the errors at la level, the strong phase 5' is set to 
be zero. 
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FIG. 2. The same as Fig.| but with 5' = vr/6 
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FIG. 3. The same as Figg but with 5' = vr/3 
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FIG. 4. The dependence of the 6 branching ratios( in units of 10~ 6 ), for 7 = 70°. The 
other parameters are at their central values. The three curves in each plot correspond to 5' = 
O(solid), 7r/6 (dashed), ir/3 (dot-dashed) The hashed bands indicate the errors(at Id) of the data. 
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FIG. 6. the branching ratio of B — > 7r°7r° (in units of 10 6 predicted as function of 5' (in degree) 
the solid line indicate the upper bound observed reported by CLEO collaboration. 
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